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Abstract 

Evolution of large-scale scalar perturbations in the presence of stiff solid (solid with pressure 
to energy density ratio > 1/3) is studied. If the solid dominated the dynamics of the universe long 
enough, the perturbations could end up suppressed by as much as several orders of magnitude. To 
avoid too steep large-angle power spectrum of CMB, radiation must have prevailed over the solid 
long enough before recombination. 


1 Introduction 

In standard cosmology, large-scale perturbations stay unchanged throughout the Friedmann ex¬ 
pansion that started after inflation, except for the last period before recombination when the 
Newtonian potential was suppressed, due to the transition from radiation to matter, by the factor 
9/10 (see, for example, [1]). The potential is not affected even by phase transitions and annihila¬ 
tions taking place in the hot universe, as long as the matter filling the universe can be regarded as 
ideal fluid. Among alternative scenarios considered in the literature there are some that relax that 
assumption, introducing a solid component of the universe formed in the early stage of Friedmann 
expansion [2, 3, 4, 5, 6, 7, 8]. The solid is supposed to have negative pressure to energy density 
ratio w, in particular, it can consist of cosmic strings (w = —1/3) or domain walls (w = —2/3). 
Such matter starts to influence the dynamics of the universe at late times only and has no effect 
on the evolution of perturbations during the hot universe period. 
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To obtain large-scale perturbations whose magnitude at recombination differs from their mag¬ 
nitude at the end of inflation, we need a solid with w > 1/3. A scenario with radiation-like solid 
(w = 1/3) was considered in [9], where it was shown that the solid produces an additional term 
in the gravitational potentials that can be large at the beginning but decays afterwards. If one 
introduces stiff solid ( w >1/3) instead, the character of the expansion of the universe changes for 
a limited period and a question arises whether this cannot cause a shift in the nondecaying part of 
the potentials, in analogy to what we observe in a universe filled with ideal fluid as it passes from 
one expansion regime to another due to a jump in w. If so, the incorporation of the solid into the 
theory, with the value of its shear modulus left free, would enlarge the interval of admissible values 
of the primordial potential, extending in such a way the parameter space of inflationary scenarios. 

A possible realization of stiff solid would be a system of equally charged particles with anisotropic 
short-range interaction. By using Yukawa potential, one obtains stiff fluid [10, 11]; however, if the 
potential is squeezed in some direction and the particles are arranged into a lattice, the system 
acquires nonzero transversal as well as longitudinal sheer modulus with respect to that direction. 

In order that a solid, radiation-like or stiff, has an effect on large-scale perturbations, the 
solidification has to be anisotropic , producing a solid with flat internal geometry and nonzero 
shear stress. Such solidification might possibly take place in case the Friedmann expansion was 
preceded by solid inflation, driven by a solid with w < 0 rather than by a scalar field [12, 13, 14, 15]. 

In the paper we study how a stiff solid formed during Friedmann expansion would influence 
the evolution of large-scale perturbations. In section 2 we derive solution for such perturbations in 
a one-component universe and establish matching conditions in a universe whose matter content 
has changed abruptly; in section 3 we determine the behavior of perturbations after the solid has 
been formed and find both nondecaying and decaying part of Newtonian potential after radiation 

prevailed again; and in section 4 we discuss the results. Signature of the metric tensor is (-1 -) 

and a system of units is used in which c = 167 tG = 1 . 


2 Perturbations in the presence of solid 

2.1 E vo hit ion equations 

Consider a flat FRWL universe filled with an elastic medium, fluid or solid, with energy density p 
and pressure p, and denote the conformal time by p and the scale parameter by a. Expansion of 
the universe is described by the equations 

r\ A 1 / 2 

a' = (g P a J , p' = -3 Hp+, (1) 

where the prime denotes differentiation with respect to 77 , 'H = a! j a and p + = p + p. 
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In a perturbed universe, spacetime metric and stress-energy tensor acquire small space-depen¬ 
dent corrections 5gp, v and ST We will use the proper-time gauge in which 5goo = 0 (the 
cosmological time t = J adg coincides with the proper time of local observers). The metric in this 
gauge is 

ds 2 = a 2 [dg 2 + 2Bidgdx l — (Sij — 2ipSij — 2E^j)dx l dx^], (2) 

where the effective equality indicates that only the scalar part of the quantity in question is given. 
Suppose the matter filling the universe has Euclidean internal geometry and contains no entropy 
perturbations. The perturbation to T is then given solely by the perturbation to g liv and the 
shift vector of matter We will use the remaining gauge freedom to impose the condition £ = 0, 
so that our gauge will be also comoving. In this gauge, the perturbation of mass density Sp = STq°, 
the energy flux density S l = — Tj° and the perturbation of stress tensor St u = ST'J are [2] 

Sp = p+{3tp + £), S' 1 = - p+B^i, St 1 -’ = - K(3ip + £)Sij - 2pE^ j . (3) 

where K is the compressional modulus, p is the shear modulus and the index ‘T’ denotes the 
traceless part of the matrix. (Our K is 2 times greater and our p is 4 times greater than K 
and p in [2]. We have defined them so in order to be consistent with the standard definitions in 
Newtonian elasticity.) 

The proper-time gauge is not defined uniquely since one can shift the cosmological time by an 
arbitrary function <5f(x). Under such shift, E stays unaltered and B and rp transform as 

B B + Sg, ip ip — HSp, 

where Sr] = a~ 1 St. This suggests that we represent B and %p as 

B = B + X , ip = ~H X , (4) 

where B stays unaltered by the time shift and x transforms as % —> % + Sg. 

We will restrict ourselves to perturbations of the form of plane waves with the wave vector 
k, B and £ oc e* k x . The action of the Laplacian then reduces to the multiplication by — fc 2 ; in 
particular, the definition of £ becomes £ = —k 2 E. For simplicity, we will suppress the factor e ik x 
in B and £, as well as in other functions describing the perturbation. They will be regarded as 
functions of g only. 

Evolution of scalar perturbations is governed by two differential equations of first order for the 
functions B and £, coming from equations = 0 and 2Goo = Too- The equations are [16] 

B’ = (3c| 0 + a- l)HB + c%\\£, £’ = -(k 2 + 2,aH 2 )B - aH£, (5) 

where a = p + /(2'H) 2 = (3/2 )p + /p, cso is the “fluid” sound speed (sound speed of the solid with 
suppressed contribution of shear modulus), c| 0 = K/p + , and cg|| is the longitudinal sound speed, 
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4 ii — 4o + (4/3 )p/p+. The only place where the shear modulus enters equations (5) is the term 
Cg||£ in the equation for B. 

Consider a one-component universe filled with a solid that has both p and p proportional to 
p. The quantity K is then proportional to p, too, since K = p+Cg 0 and c | 0 = dp/dp. Mechanical 
properties of such solid are given completely by two dimensionless constants w = p/p and p = p/p. 
To simplify formulas, we will often use the constant /3 = p/p+ = p/w + , where tu + = 1 + w, instead 
of p. 

For constant w and p, the quantities appearing in the equations for B and £ are all constant, 
except for the Hubble parameter that is proportional to p~ x ■ Explicitly, 

3 4 

a = 2 W +’ 4o =u h c 2 s\\= w + -I3 = w, H = 2up~ 1 , 

where u = 1/(1 + 3 w). With these expressions, equations for B and £ simplify to 

B' = u(l + 9w)p~ 1 B + w £, £'=—{k 2 -\-lS,u 2 w + p~ 2 )B — Zuw + p~ l £, ( 6 ) 

and after excluding £, we arrive at an equation of second order for B , 

B" + 2vp~ 1 B' + [q 2 - (2v - b)p~ 2 ]B = 0, (7) 

where q = \/duk, v = u{ 1 — 3 w) and b = 2Au 2 p. The equation is solved by Bessel functions 
of the argument qp, multiplied by a certain power of p. We are interested only in large-scale 
perturbations, that is, perturbations stretched far beyond the sound horizon. Such perturbations 
have qp <C 1, hence we can skip the term q 2 in the square brackets in (7) to obtain 

B = p(cjp~ m + c Y p~ M ), (8) 

where the parameters m and M are defined in terms of the parameters v = v + 1/2 = (3/2)w(l — w) 
and n = V v 2 — b as m = v — n and M = v + n. The constants are denoted cj and cy to remind 
us that the two terms in ( 8 ) come from the Bessel functions J and Y. 

The function B is non-oscillating for b < v 2 and oscillating for b > v 2 . Solutions of the second 
kind are well defined if the solid was not present in the universe from the beginning, but was 
formed at a finite time. Here we will restrict ourselves to the solutions of the first kind, which 
means that we will consider only values of the dimensionless shear stress p < (3/32)(l — w) 2 . 

An approximate expression for £ is obtained by inserting the approximate expression for B 
into the first equation in ( 6 ). In this way we find 

£ = Cjp~ m + Cyp~ M , (9) 

where cj and cy are defined in terms of cj and cy as cj = — (l/th)(3/2 — n)cj and cy = 
— (l/u>)(3/2 + n)cy. 
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2.2 Potentials $ and \[r 


Scalar perturbations we are interested in are most easily interpreted in the Newtonian gauge, in 
which the metric is 

ds 2 = a 2 [(l + 2^)drf - (1 - 2 'F)dx 2 ]. ( 10 ) 

Let us express the potentials $ and 4' in terms of the functions B and £. If we perform explicitly 
the coordinate transformation from the proper-time to Newtonian gauge, we find (see equation 
(7.19) in [1]) 

’F = H{B- E'). (11) 

For d> we could proceed analogically, but it is simpler to use Einstein equations. If we write the 
scalar part of the stress tensor as a sum of pure trace and traceless part, r' <J = r^Stj + t- 2 ^, 
from equations 2 Gij = T, l3 we obtain that the difference of d> and dt is given by the latter quantity 
(see equation (7.40) in [1]), 

A$ = d> - d> = -r^a 2 . 

2 

By inserting here from the third equation (3) we find 

Ad> = —fia 2 E. (12) 

We can see that in a universe filled with an ideal fluid (/j = 0) the potentials $ and dt coincide. 

After inserting into the expression for d' from the second equation in (5) and into the expression 
for Ad> from the first equation in ( 1 ), we arrive at 

<F = ~k~ 2 aU 2 {mB + £), A$ = 6 Jik~ 2 H 2 £. (13) 

For the one-component universe introduced before, expressions for dt and Ad> become 

= — 6u 2 w + (kr])~ 2 (6ur]~ 1 B + £), Ad> = 24 u 2 £. (14) 

With B and £ given in (8) and (9), both $ and dt are linear combinations of g^ 2 ^ m and rj~ 2 ~ M , 
For an ideal fluid m = 0 and M = 2v, so that we expect the function $ to be linear combination 
of g~ 2 and ri~ 2v +, where u + = 1 + u. This is, however, not true because the coefficient in front of 
rj~ 2 turns out to be zero. Thus, if we want to establish how >F looks like for an ideal fluid, or how 
d> and d> look like for a solid with small jl, we must add the next-to-leading term to the J-part of 
both expressions (8) and (9). The term is suppressed by the factor ( qr]) 2 , therefore the J-part of 
d> for an ideal fluid is constant and the J-part of $ and d> for a solid with small jl acquires a term 
proportional to r/~ m . For a universe filled with an ideal fluid we have 

B = d{cj + c Y r]^ 2v ), £ = cj + c Y r]^ 2v , (15) 
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where cj and cy are defined in terms of cj and cy as cj = —6ucj and cy = —3 u(wy/w)cy. After 
computing the additional terms in B and £ and inserting the resulting expressions into equations 
(14), we arrive at 

< f ) = Cj + Cyr]~ 2v+ , (16) 

where Cj and Cy are defined in terms of cj and cy as Cj = 3 u 2 (w+/v+)cj and Cy = 
12 u 2 w + vq~ 2 cy. 


2.3 Transitions with jump in w and ft 


Suppose the functions w v and p v change at the given moment p t r (“transition time”) from (wi, pi) 
to (w 11 , fin) = (wi + A w,pi + A p). (We have attached the index rj to the symbols w and p in 
order to distinguish the functions denoted by them from the values these functions assume in a 
particular era.) Rewrite the first equation in (5) as 


B' = - 2 


§ 0 (3 TLB + £) + (^—w v + — l^TLB + 


where 


2 dp dw v 

c s o = T~ = w v + P- 


(17) 


(18) 


dp dp 

Because of the jump in w v , there appears (5-function in c| 0 , and to account for it we must assume 
that B has a jump, too. However, on the right hand side of equation (17) we then obtain an 
expression of the form “0-function x (5-function”; and if we rewrite B' as 

dB , d.B 

B = —p = -3 npw v+ —, 

dp dp 

on the left hand side there appears another such expression. To give meaning to the equation we 
must suppose that w v changes from w / to wu within an interval of the length Ap <C p tr , and 
send A p to zero in the end. If we retain just the leading terms in equation (17) in the interval 
with variable w, we obtain 

dw r i 


dB _ _ / t tr \ ClWr, 
do V 3 T-Ltr-J do 


W " + Tp = 


4 

iditr' 


dp 


(19) 


where we have used the fact that, as seen from the second equation in (5), the function £ is 
continuous at r) = rj tr . The solution is 

3rLtr U7 

Denote the jump of the function at the moment r] s by square brackets. To determine [£>], we 
express Bi and Bn in terms of wi+ and wii+, compute the difference Bn — Bi and use the 
expression for £>/ to exclude C. In this way we find 


wu .l V 37t+ r / 




( 20 ) 
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Note that the same formula is obtained if we assume that the functions with jump are equal to 
the mean of their limits from the left and from the right at the point where the jump occurs. 

To justify the expression for [£>], let us compute the jump in U/. It holds 

[*] - — ^k~ 2/ H 2 r (SHtr[w v+ B] + A w£tr), 

and if we write [w^+B] = vj n+[B] + A wBi and insert for [B\, we find that [\&] vanishes. This must 
be so because for T we have (see equation (7.40) in [1]) 

T" + H{ 2'L' + $') + (2 n' + H 2 )^ = -ijrW, 

where the bar indicates that the quantity is computed in Newtonian gauge. A jump in T 
would produce a derivative of (5-function in the first term, but no such expression with opposite 
sign appears in the other terms. 

The jump in B’ can be found from equation (17) by computing the jump of the right hand 
side, with no need for the limiting procedure we have used when determining the jump in B. The 
result is 

[B’] = 4 —UtrBtr + ( 5 ~ SWlI Aw + | A p)s tr . ( 21 ) 

wu+ V o wn+ 6 J 

3 Scenario with stiff solid 

3.1 Expansion of the universe 

Suppose at some moment r\ s the hot universe underwent a phase transition during which a part 
of radiation (w = 1/3) instantaneously turned into a stiff solid (w > 1/3). In a one-component 
universe with given parameter w, the density of matter falls down the faster the greater the value 
of vj. As a result, if the solid acquired a substantial part of the energy of radiation at the moment 
it was formed, it dominated the evolution of the universe for a limited period until radiation took 
over again. Let us determine the function a(rf) for such universe. 

Denote the part of the total energy that remained stored in radiation after the moment ?y s by 
e. In the period with pure radiation (?; < ij s ) the mass density was p = p s (a s /a) 4 , so that from 
the first equation in ( 1 ) we obtain 

/I N 1 / 2 

« = C V , C=[- Ps a 4 s ) . (22) 

In the period with a mix of radiation and solid (77 > r] s ) the mass density is 

p = ep s (a s / a ) 4 + (1 - e)p s (a s /a) 3w+ = p s (a s /a) 4 [e + (1 - e)(a s /a) A ], 
where A = 3w + — 4. As a result, the first equation in (1) transforms into 

a' = C'[e + (l-e)(a s /a) A ] 1/2 . (23) 
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For w > 1/3 the parameter A is positive, therefore the second term eventually becomes less than 
the first term even if e <C 1 . 

Suppose radiation retained less than one half of the total energy at the moment of radiation-to- 
solid transition (e < 1/2). The subsequent expansion of the universe can be divided into two eras, 
solid-dominated and radiation-dominated, separated by the time i^ ra d at which the mass densities 
of the solid and radiation were the same. The value of rj ra d is given by 

a-rad = a s (e -1 - 1) 1/A . (24) 

Suppose now that the post-transitional share of energy stored in radiation was small (e <C 1). 
The universe then expands by a large factor between the times rj s and r] ra d, 

drad = d s C ^ tt s , 


and can be described in a good approximation as if it was filled first with pure solid and then with 
pure radiation. Thus, equation (23) can be replaced by 


f C(a s /a) A/2 for 77 < r] ra d 
l s/eC for 77 > ?7 

rad 


(25) 


The solution is 

. f [(A/2 + 1 )af /2 Crj] A/2+1 for 77 < r] rad 
a = < 

1 \ftCf) for 77 > 1J r ad 

where fj and 77 are shifted time variables, fj = 77 — 77 * and fj = 77 — 77 **. 
expression for a ra d we obtain 


(26) 


From the approximate 


1 


?7 rad — 


A/2 + 1 

and by matching the solutions at rj s and 77 ra d we find 


A / 2+1 

e A Vs 


A, 


A/2 

V* ~ /\J2 + ;L^ S ’ ? ?** — T^Vradj 


(27) 


(28) 


Note that equation (23) solves analytically for w = 2/3 and w = 1, when A = 1 and A = 2. 
We do not give these solutions here since will not need them in what follows. 

The two equations in (28) can be rewritten to formulas for the ratios of shifted and unshifted 
times, 


Vs _ 1 _ U fjrad _ A | ^ _ u 0 

7? s A/2 + 1 U 0 ' fjrad 2 U ’ 

where uq is the value of u in the radiation-dominated era, uq = 1/2. These equations stay valid 
also after we replace radiation by an ideal fluid with an arbitrary pressure to energy density ratio 
wq. To demonstrate that, let us derive them from the condition of continuity of the Hubble 
parameter. If the universe is filled in the given period with matter with the given value of w, its 
scale parameter depends on a suitably shifted time fj as a oc fj 2u . Thus, its Hubble parameter is 
T~L = 2ufj~ 1 and the requirement that TL is continuous at the moment when w changes from wi to 
wu is equivalent to fjn/fji = uu/ui. 



3.2 Behavior of the function B 


We are interested in large-scale perturbations in a universe in which the parameters w and p 
assume values (wq, 0 ) before rj s , ( w,p ) between p s and r) T ad, and (wq, 0) after rirad . (Most of the 
time we will leave wo free, only at the end we will put Wq = 1/3.) Denote the functions describing 
the perturbation before p s and after rj ra d by the indices 0 and 1 respectively, and keep the functions 
referring to the interval between and rj ra d without index. If only the nondecaying part of the 
perturbation (the part with constant <h) survives at the moment p s , the functions Bo and So can 
be replaced by their J-parts, 


Bo = cjqt], Sq = cj o = — 6 mqCjo- 


(29) 


For the functions B and S we have expressions ( 8 ) and (9) with ?/ replaced by rj and for the function 
B\ we have the first equation (15) with cj and c Y replaced by cji and cyi, v replaced by vq and 
77 replaced by rj. All we need to obtain the complete description of the perturbation is to match 
the expressions for Bo, B and B\ with the help of the expressions for So and S at the moments rj s 
and ri ra d■ 

At the moment r] s , the jumps in w v and are A w s = w — wo = Aw and Ap s = p. By using 
these values and the identity So = —SHsBos, we find 

[B\s = 0 , [B'] s = -(^Aw-^)s 0 , 

Denote Xq = cjq- Equations for the unknowns x = CjrjJ m and y = Cyfj~ M are 


x + y=— xq, (1 - m)x + (1 - M)y = 

u 


1 + 


8 «°(^Aw - (3^j 


Xo, 


(30) 


and their solution is 


x = — ^-(M - 8uf3)xo, y = - — - 8uf3)x 0 - 

u 2 n u 2 n 


(31) 


At the moment r] ra d, the jumps in w v and fi r) are A w ra d = — A w and A (3 ra d = —p,. By inserting 
these values into the expressions for [B] and [ B 1 ] we obtain 


[B]rad = — (Brad + , [B']rad = -4 — UradPrad ~ ( 

Wox V 3 Tirad'' Wox V 


Aw 
W 0 + 

Introduce the constants 


Aw , 
W0+ 


' 5 — 3wq 
6 w 0 + 


Aw ■ 


■'rad’ 


X = cjrj; a m d = P~ m X, Y = c Y fj^ = p- M y, (32) 

where p is the ratio of final and initial moments of the period during which the solid affects the 
dynamics of the universe, p = fjrad/ijs■ Equations for the unknowns x = cj\ and y = cyifj~^° are 

x + y= — (KjX + KyY), i +(1-2^ )y = LjX + L Y Y, (33) 

u 0 
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where the coefficients on the right hand side are defined as 


Kj 


W 0+ L 


W+ 


Aw 
6 uw 


(in + 6 uw) , 


Ky = ditto with m —> M, 


and 


8uAw m + 6uw/'5 — 3wo 4 \ 

Lj = l — m -1-1- Aw H— jj, , Ly = ditto with m —> M, 

ieo+ w V 6u>o+ 3 / 


6 w 0 + 


The solution is 


~x= ^-(MjX + M y Y). § = -J-(NjX + NyY) 
2v 0 


(34) 


with the constants M a and N a , a = J, Y, defined in terms of the constants L a and K a as 


u u 

M a = L a - (1 - 2u 0 ) — K a , N a = L a - -K a . 

U 0 UO 


3.3 Behavior of potentials 

Knowing how the function B looks like, we can establish the time dependence of the Newtonian 
potential <k and the potential describing the curvature of 3-space \1/. Before the time rj s , both 
potentials are the same, 4>o as well as tko = Cjq ~ xq. Between the times ?y s and i] ra d, the 
potentials are given by the two equations in (14) with r/ replaced by f). With B and £ inserted 
from equations (8) and (9), both $ and tk become sums of terms proportional to and 

fj~ 2 ~ M . We have already mentioned that for fi = 0 the coefficient in the first term in 4> = tk is 
zero, and one easily verifies that for w > 1/3 and fi close to zero the first coefficient in both $ and 
\k is proportional to /2. (After a simple algebra we find that it is proportional to m(l — 4/3) — 8/3 
and m — 8w/3 for $ and \k respectively, with m reducing to b/(2v) = 8uw+f3/(l — w) in the limit 
/3 <C 1.) The coefficients contain the constants cj and cy and if we use cj (X x and cy oc y with 
x and y given in equation (31), we find that the second coefficient is proportional to fi, too. (In 
the expression for y we encounter the factor m — 8iif3 again.) Both coefficients contain also the 
factor xo ~ <l>o, therefore for 77 close to rj s we have <k as well as tk ~ /i(fc?)) _2 4> 0 . As 77 grows, the 
first correction to the term proportional to fj^ 2 ~ m , which is of order 4>o, may take over while the 
perturbation still remains stretched over the horizon. However, in order that our approximation is 
valid, this term must be negligible in the first period after the moment rj s . (Note that this does not 
hold for the potential \k just after r) s : it equals 4>o at r/ s , hence it is dominated by the correction 
term for a short period afterwards.) As a result, /2 must be not too close to zero, (1 ( kfj s ) 2 . 

For large enough fi, $ and T can become much greater in absolute value not only than $ 0 , 
but also than 1. The theory then seems to collapse, but it does not because, as can be checked 
by direct computation, kB, ip and £ remain much less than 1. (A detailed discussion for Aw = 
0 can be found in [17].) Thus, the proper-time comoving gauge which we have implemented 
instead of more common, and intuitively more appealing, Newtonian gauge, is not only convenient 


10 



computationally, but also preferable on principal grounds. Without it we would not know that the 
perturbations stay small and the linearized theory stays applicable after a solid with above-critical 
parameter ft, was formed, causing the potentials $ and T to rise beyond control. 

We are interested in the potential $ after the moment r) ra d, when both potentials coincide 
again. Denote the nondecaying part of $ in that period as It holds $1 , nd = Cji, and by 

using the relation between Cj and cj we obtain 


* _ O 2 w 0+ ~ 

Ql,nd — oUq X. 

^0+ 


(35) 


Here we must insert for x from equation (34), with X and Y given in equation (32), x and y given 
in equation (31) and Xq given by 


q 2 w o+ 

$0 = 3u 0 - Xq. 

Vq+ 


The resulting expression for <nd is 


<t>l,nd = ^--^( Mjp~ m - M Y p- M )$ 0 , 

2 i/q u 2n 


(36) 


with the coefficients Mj and My defined as 


Mj = Mj (.M - 8 u/3 ), M y = M Y (m-8u/3). 


After some algebra the coefficients reduce to 


Mj = 2uq — M — 6 , My = ditto with M 
uq 


TO. 


(37) 


Let us now determine how fast the function $ approaches its limit value. Denote z = qofj, 
where qo = Wok. The decaying part of 4> in the period under consideration is 


A4>! = -2vo + 1 ^±-(Njp~ m - N Y p~ M )~z;* d C- 2 " 0 +$o, 


(38) 


where £ is rescaled time normalized to 1 at the moment rj ra d , C = l/fjrad, and the coefficients Nj 
and N y are defined in terms of Nj and Ny in the same way as the coefficients Mj and My in 
terms of Mj and M Y . After rewriting the former coefficients similarly as we did with the latter 
ones, we obtain 


Nj =N Y = - 


w 0 
^0+ 


26. 


(39) 


From these equations and equations (36) and (37) we find that the ratio of the decaying and 
nondecaying part of $ at the moment of solid-to-radiation transition is 


A$i 


$ 


1 ,nd 


= R5 r ad’ R = 4 ^0 D>+ W ° 


2uob 


rad 


wo + 2uQu[ncoth(nlogp) + v] — uob 


(40) 


The ratio is greater than one for (3 > z^ ad . The function $ is then dominated by the decaying 
term at the moment p ra d, the nondecaying term taking over later, at the moment rj n d given by 

. 1 — 3 - 


z nd = R^o+ z rad 


(41) 
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The exponent at z ra d is positive for any wo < 1 (it equals 1/3 for Wq = 1/3) and the constant R 
is of order 1 or less. Thus, if the perturbation was stretched over the horizon at the moment the 
fluid originally filling the universe started to be dominating again ( z ra d -C 1), it will stay so at the 
moment the nondecaying term prevails over the decaying one (z„d 1). 

The time ri ra d must not be too close to the time of recombination r) re , if the spectrum of 
large-angle CMB anisotropies is not to be tilted too much. If we denote the wave number of 
perturbations crossing the sound horizon at recombination as k^°\ the perturbations with the 
longest wavelength that can be observed in CMB have k ~ O.OIP 0 / For w o = 1/3, Newtonian 
potential after the moment r] ra d is $i = (1 + Rz~ ad C, s )^\ tn d = (1 + Rz ra dZ~ 3 )^i tn d, and if we take 
into account that the value of z$) is approximately 1, we find 


(fc(°)) 

r = 


$i(0.01fc(°)) 


1 + Rz, 


(0) 

rad 


i + 


_= i _ in 4 r|(°) 


The observational value of r is 0.01" s_1 , where ng is the scalar spectral index , a characteristic of 
perturbations whose deviation from 1 (about —0,04 according to observations) describes the tilt 
of the scalar spectrum. If we allow for a tilt of the primordial spectrum, too, the expression for 


r must be multiplied by 0.01 nso 


Denote p* = 1 /zf^ d = rjre/rjrad = 


e/arad — Trad/Tre and 


require that ns differs from nso at most by some A ns <C 1. To ensure that, p* must satisfy 


p* > 2 x lO^iJArig 1 . 


(42) 


For numerical calculations we need the value of p. It is a ratio of times, but can be rewritten 
in terms of a ratio of scale parameters or temperatures, P = a ra d/a s = T s /T ra d, as 


P = P 2" 


(43) 


The value of p, or equivalently, P, determines the interval of admissible w’s. To obtain it, note 
that for wq = 1/3 equation (24) yields P = (e -1 — l) 1 ^ = e -1 / A , or 


P = e si™. 


(44) 


(This is consistent with equation (27), which can be rewritten as p = e / a + = e 6 «a«.) Thus, 
the jump in the parameter w for the given ratio P must satisfy 


Aw = 


log 1/e 
3 logP 


1 

3 log P ’ 


(45) 


The dependence of the quantities <j> = ^i, n d/^o and R on the parameter /3 is depicted in fig. 1. 
The values of w o and w are 1/3 and 2/3 on both panels and the solid and dotted lines correspond 
to P = 10 3 and P = 10 13 respectively. The lines are terminated at 0 = 1/160, which is the 
maximum value of /3 admitting non-oscillating solutions in a solid with w = 2/3. 
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Fig. 1: Final value of Newtonian potential in a universe with stiff solid (left) and normal¬ 
ized ratio of decaying to nondecaying part of the potential at solid-to-radiation 
transition (right), plotted as functions of dimensionless shear modulus 


For completeness, we have included also values /3 < 0 into the graphs. The transversal sound 
speed squared is negative for such /3, so that the vector perturbations start to grow exponentially 
once they have appeared. As a result, the theory is acceptable only if such perturbations are 
produced neither during inflation (which is the case in simplest models) nor in the subsequent 
phase transitions. 

The parameter P assumes the smaller value if, for example, the solid dominated the dynamics 
of the universe between the electroweak and confinement scale, and the greater value, if the solid 
was formed as soon as at the GUT scale and dominated the dynamics of the universe up to the 
electroweak scale. Unless the parameter w of the solid is close to that of radiation, the fraction of 
energy which remains stored in radiation after the solid has been formed must be quite small in 
the former case and very small in the latter case. For w = 2/3 this fraction equals 1/P, so that 
for the greater P the mechanism of the radiation-to-solid transition must transfer to the solid all 
but one part in 10 trillions of the energy of radiation. 

The quantity cj) is the factor by which the value of the potential $ changes due to the presence 
of stiff solid in the early universe. From the left panel of fig. 1 we can see that $ is shifted upwards 
for /3 < 0 and downwards for (3 > 0, and the enhancement factor decreases monotonically with 
/3, the steeper the larger the value of P. For maximum /3 the function $ is suppressed by the 
factor 0.41 if P = 10 3 and by the factor 0.004 if P = 10 13 . The quantity R determines, together 
with the parameter A ns, the minimal duration of the period between the moment when radiation 
took over again and recombination. According to the right panel of fig. 1, the temperature at the 
beginning of this period had to be at least 8 x 10 3 Ang 1 P re = 0.2 (Ans/0.01) -1 MeV for maximum 
/3 and P = 10 3 . 


13 






4 Conclusion 


We have studied a scenario with stiff solid appearing in the hot universe and dominating the 
evolution of the universe during a limited period before recombination. In comparison with the 
scenario containing radiation-like solid [9], a new effect is that the nondecaying part of Newtonian 
potential becomes suppressed. This might raise hope that the tensor-to-scalar ratio is enhanced, 
which would surely be interesting from the observational point of view. However, a straightforward 
calculation shows that the tensor perturbations are suppressed by exactly the same factor as the 
scalar ones. The shift in Newtonian potential towards less values means that the rms of primordial 
potential was in fact greater than supposed. As a result, there appears an additional freedom in 
the choice of the parameters of inflaton potential; for example, one can use potentials with smaller 
inclination of the plateau than in the case without solid when implementing slow-roll inflation. 
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